Computational Understanding and Manipulation of Symmetries 
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Abstract 

Symmetry structure in natural and artificial systems 
such as crystallography, chemistry, physics, and permu- 
tation puzzles, et cetera, can facilitate understanding and 
manipulation of these systems, as is well-known in the 
mathematical sciences and from the algebraic theory of 
groups (the latter now with about two centuries of con- 
tinuous development). However, until now, these meth- 
ods have not been fully exploited in Artificial Intelli- 
gence (AI). 

We show how AI systems based on computer alge- 
braic methods can make use of such mathematical sym- 
metry structure to automatically generate and manipu- 
late hierarchical coordinate systems for finite systems 
whose generating symmetries are given. Such hierarchi- 
cal coordinate systems correspond to subgroup chains 
in the group structure determined by generating sym- 
metries of the system. We demonstrate how, for any 
finite-state symmetry system, these coordinate systems 
(1) can be generated automatically, i.e. deriving formal 
models for understanding the finite-state symmetry sys- 
tem, and (2) can be deployed in manipulating the system 
automatically. Thus, without learning, manipulation of 
such systems is reduced via algebra to sequential com- 
putation in a hierarchy of simple (or simpler) group co- 
ordinates. 

This general, implemented method is illustrated with 
examples from coordinate systems on permutation puz- 
zles such as Rubik's cube. 



1 Introduction 

Algebra and the theory of permutation groups is 
well-known from applications of groups in chem- 
istry and crystallography [Raman, 2004], in physics 
| Iones, 1990] , and more recently of semigroups 
and groups in systems biology, genetic regulatory 



networks, and biochemistry [Eg ri-Nagy et al, 2008 
Rhodes, 20 09]. Unlike machine learning or optimiza- 
tion techniques, algebraic machine intelligence can 
without any learning derive coordinate systems not only 
on states of a structure, but also on its transformations, 
i.e. on operations for manipulating the structure. These 
are unlike the methods, spatial structures, or cogni- 
tive and semantic maps presently exploited for theorem 
proving, path-planning, automated reasoning, but simi- 
lar to coordinate system occurring elsewhere in science. 
Coordinate systems, like the ones studied here on gen- 
eral symmetry structures from the viewpoint of AI ap- 
plications using computer algebra, also arise in conser- 
vation laws in physics. As Emmy Noether showed in 
the first part of the 20th century, invariants preserved 
under conservation laws correspond exactly to group- 
theoretic symmetry structures in physics [Olver, 2 000 1; 
moreover, such invariants for physical systems give rise 
to coordinate systems of exactly the type described here 
| Nehaniv, 19961 |Rhodes, 20091 . 



Algebra, Cognitive Modeling, Coordinate Systems 

AI techniques for cognitive modeling, machine learn- 
ing and optimization so far have made relatively lit- 
tle use of abstract algebra. Cognitive architectures 
such as SOAR, ACT-R and SAC have been applied 
to build AI systems that model human cognitive ca- 
pacities more or less intending to emulate faithfully 
the structure of cognition in humans, with applica- 
tions ranging from autonomous control of aircraft based 
on subsymbolic rule-extraction, to human-style learn- 
ing of arithmetic or natural language, to predictive 



evaluation of user interfaces (e.g. [Laird et al., 19871 , 
I IReder and Schunn, 1996| , QLi et al, 20031 ). But as- 
pects of human and machine cognition can also involve 
understanding of hierarchical processes with dynami- 
cal structure as evidenced by the object-oriented meth- 



ods [Nehaniv, 1997] or the place-value representation 



in human number systems [Zhang and Norman, 1995], 
which have close connections to algebra. Here we are 
interested in the study of AI models that can derive, 
represent, and manipulate this type of knowledge, but 
without necessarily seeking to model human capaci- 
ties faithfully. Models for understanding of finite-state 
(and more general) dynamical systems phenomena in 
general exhibit such feedforward, coarse-to-fine, hierar- 
chical structure related to algebraic coordinate systems 



metry structures (called '[permutation] groups' 
in algebra) via the Frobenius-Lagrange embed- 
ding. Related work |Egri-Nagy and Nehaniv, 2004 



I Krohn et al.,T968l|Nehaniv, 1997[|Rhodes, 20091 . Our 
work shows how such coordinate systems can be derived 
and exploited automatically. 

By a coordinate system on a symmetry structure (or 
more general structure), we mean a notational system in 
the broadest possible sense, with which a human or arti- 
ficial agent can address building blocks of the structure 
and their relations in a decomposition, thus gaining a 
convenient way for grasping the structure of the original 
phenomenon and possibly getting tools for manipulating 
the components. 

An obvious example is the Descartes coordinate sys- 
tem, where we can uniquely specify any point of the 
n-dimensional space by n coordinates. However, this is 
an example of a spatial and inherently non-hierarchical 
coordinate system for a totally homogeneous space. In 
general, different coordinates have different roles, ad- 
dressing 'parts' of the system different in size, function, 
etc. The natural example of a hierarchical coordinate 
system is our decimal positional number notation sys- 
tem: different coordinates correspond to different mag- 
nitudes. The examples also show that a coordinate sys- 
tem is very much the same thing as a decomposition: 
the space is decomposed into dimensions, an integer is 
decomposed into ones, tens, hundreds, etc. 

From this viewpoint coordinate systems become cog- 
nitive models, means of knowledge representation. The 
above examples show the usefulness of these coor- 
dinates, but how can we obtain these models? The 
good news, and the main promise of this research di- 
rection, is that we can get them automatically! In 
algebraic automata theory, the Prime Decomposition 
Theorem says that every finite state automata can 
be decomposed into a hierarchical coordinate system 



DKrohnetaL, 1968|[EiTehberg, 1976||Holcombe, 19821 . 
Therefore the way of representing knowledge becomes 
algebraic, semigroup- and group-theoretical, which is 
really very different from other, well-established AI 



methods mainly based on logic (e.g. [Ne well, 1990 



Egri-Nagy and Nehaniv, 2008] |Egri-Nagy et al, 2008) 



has computationally implemented the automated gen- 
eration of transformation semigroup decompositions 
along the lines of the Krohn-Rhodes theorem, but did 
not pursue the 'simpler' problem of obtaining coordi- 
nate systems on permutation groups as we do here, and 
hence our work is both complementary and necessary 
for providing a full decomposition in the more general 
setting where operations need not be invertible. 

The use of algebra in cognitive modeling also occurs 
in other non-traditional applications. For example, the 
approach by Fauconnier-Turner-Goguen to conceptual 
blending and metaphor [F auconnier and Turner, 2 003], 
uses category-theoretic pushout computations with 
computer algebraic implementations by Goguen to au- 
tomatically generate conceptual blends and metaphors 
relating two (or more) conceptual domains (knowledge 
of which is modeled by small categories), as well as 
applications to the semiotics of user-interface design, 
and formal specification for imperative programs built 
up and verified in a hierarchical manner [Go guen, 1999| 
Goguen and Malcolm, 2000 1. 

As there are many different ways of understanding of 
the same thing, there are many different coordinate sys- 
tems for the same structure. Some of them may be in- 
tuitive for humans, while others will suit computational 
manipulation better, so the range of intelligent 'users' of 
the coordinate system is not restricted. 

Although much of the mathematical theory required 
here is very old, the proper computational tools were 
missing, therefore the idea of hierarchical coordinate 
systems giving understanding computationally has not 
been much studied nor applied. Now such tools are 



available open-source | Egri-Nagy and Nehaniv, 2004 



Brachman and Levesque, 20 04 ] ) . 

Here we concentrate on coordinatizing sym- 



Egri-Nagy and Nehaniv, 200SJ for transformation semi- 
groups, complemented by our work reported here for 
groups giving fine detail on coordinatizing the groups 
involved. The mathematical significance of these coor- 
dinate systems is immediate, but as they capture one of 
the basic aspects of our cognitive capabilities, namely, 
hierarchical representation, they might also play signif- 
icant role in Artificial Intelligence. 

2 Mathematical Preliminaries 

A function p : X — > X on the set X is called a permuta- 
tion if it is one-to-one and onto (invertible). A permuta- 



Hon group (X, G) is a set G of permutations closed un- 
der composition (multiplication, usually denoted by ■), 
together with the state set X on which the mappings act. 
It is called a symmetry group if certain structure on X 
is preserved by all p G G. If it is clear form the context 
we omit the state set and write simply G. For x G X 
and g G G, we write x ■ g for the result of applying 
permutation g to x. The action is faithful if whenever 
x ■ g = x ■ g' holds for all x G X then 5 = </. The 
group contains the identity 1 and inverse map for 
each element p, thus everything can undone within a 
group. A group G acting on X by permutations need 
not necessarily be faithful, i.e. need not be a permuta- 
tion group. The group consisting of all permutations of 
n points is called the symmetric group S n , while C n de- 
notes the cyclic group on n objects permuted cyclically. 
A subset of a group called a subgroup if it is closed un- 
der inversion and the group's multiplication. If (X, G) 
is faithful, it is a naturally a subgroup of the symmet- 
ric group on X. For more on elementary group the- 
ory see for instance QHall, 19591 IRobinson, 1995| , and 



on permutation groups see [Dixon and Mortimer, 1996 



Cameron, 19991 . As is standard, we use cycle notation 
to denote permutation group elements, e.g. (1 3)(2 4 5) 
denotes the permutation swapping 1 and 3, cyclically 
taking 2 to 4, 4 to 5, and 5 to 2, while leaving any other 
objects fixed. 

As a special case, a group can act on itself, i.e. the 
group elements are the states and each g$ G G maps 
G — > G by right multiplication, g 1— > g ■ go- This is 
called the right regular representation, and it enables us 
to identify the group element with its effect, which will 
be really handy for permutation puzzles. However, this 
representation poses problems computationally due to 
the possibly large size of the state set. 

Given permutation groups (X, G) and (Y, H) their 
wreath product is the permutation group 

(X, G) l (Y, H) = (X x Y,G x D) 

where D = H x is the set of all possible functions from 
X to H. A state in the wreath product is expressed by 
two coordinates (x,y), x G X,y G Y. The group ele- 
ments are coordinatized similarly by (g, d) where g G G 
and d : X — > H is the dependency function, a 'recipe' 
to find an element that should be applied on the second 
level based on the (previous) state of the first (top) level. 
Thus a permutation of X x Y is given by: 

(x,y) ■ (g,d) = (x- g,y- d(x)) 

from which we can see that on the top level the action 
is independent from the bottom level, but not the other 



way around, hence the hierarchical nature of the wreath 
product. Wreath products are generally huge structures, 
but in practice we deal with some substructures with the 
dependency functions limited. The wreath product is 
easy to generalize for more levels. 



3 Lagrange Coordinatizations of 
Groups 

The basic idea of the Lagrange Decomposition is that 
given a subgroup H of G, we form the set of cosets 
G/H = {Hg : g G G}, i.e. the subgroup H and its 
translates within G; these partition G, and G acts on 
G/H by right translations: Hg Hgg^. Moreover, 
we do not need to act on whole cosets but on arbitrary 
but fixed representatives of cosets (g G Hg = Hg). The 
action may not be faithful, so we denote G made faithful 
by 

Theorem 1 (Lagrange Decomposition) Let (X, G) be 

a transitive permutation group and (X, H) be a sub- 
group of it. Then (X, G) admits the following coordina- 
tization 

(G/H,G)l(H,H) 
corresponding to the subgroup chain G > H > (1). 

Thus given a state x G X we will coordinatize it by 
x = (x\,X2), where each x% is a coset representative 

( Xl eG/H,x 2 eH/(l)). 

By refining the underlying subgroup chain we can 
make the component groups much simpler or simple 
(i.e. having only trivial homomorphic images). This al- 
lows one to iterate the Lagrange coordinatization so the 
problem of understanding the permutation group (X, G) 
is reduced to understanding much simpler permutation 
groups linked up in a feedforward manner. Therefore 
getting a coordinatization corresponds to devising a sub- 
group chain. 

For building subgroup chains certain subgroups are 
very useful. The stabilizer G a is the subgroup of G, 
which fixes x G X. Point-wise and set-wise stabilizers 
can be defined for sets of states as well. By iterating and 
refining the Lagrange construction, we have 



'Here G is a quotient group G/K, where K is the core of H in 
G, that is, the largest normal subgroup of G contained in H. Thus 
K — P| gG g~ 1 Hg. See standard references [Robinson, 1995 
ICameron, 1999) . 



Theorem 2 (Frobenius-Lagrange Coordinatization) 

Let (X, G) be a transitive permutation group and let 
G = G\ > G2 > • • • > G n+ i = G a be a subgroup 
chain ending at the stabilizer G a of some state a € X. 
Then, in the notation above, (X, G) is coordinatized by 
embedding in the wreath product 

(G1/G2, G\) l ■ ■ ■ l (G n /G n+ i,G n ) 

where Gi is Gi modulo the core of Gj+i in Gi. 
Moreover, since (X, G) corresponds isomorphically to 
(G/G a , G), the states have coordinates given by the bi- 
jection 

G a x n • • • x\ < £v- (jcij . . . , x n ) 

where each Xi is the fixed coset representative ofGi+iXi 
in Gi. 

Remarks. (1) The number of coordinate tuples in The- 
orem 12 is exactly \X\. Each state has a unique coor- 
dinatization. (2) In Theorem Q] the number of possible 
coordinate tuples is \G\ = \X\ • \G a \, for any a € X, as 
each point has exactly | G a | different possible coordina- 
tizations. 



4 Coordinate Manipulation 

The Frobenius-Lagrange decomposition gives us, in 
terms of a coordinate system, a structured view of the 
group, i.e. we can address its parts conveniently and 
with arbitrary precision. However, we would like to use 
the coordinate system dynamically, not just the statical 
view. We would like to calculate with it, finding ma- 
nipulative operations taking one state to another desired 
state, or, equivalently, from one tuple of coordinate val- 
ues to another one using the elementary symmetry oper- 
ations of the original structure. 

4.1 Component Actions 

For establishing the connection between the original 
group and the coordinatized one we need to have a way 
to express a permutation as coordinate actions. Given 
a group element g G G and a coordinatized state x = 
(xi, . . . , x n ), we can calculate the coordinatewise com- 
ponent actions by the following recursive calculations: 



91 ■= 9 



gi is not the identity. Note the hierarchical structure: gi 
depends only on g and (x±, . . . , 2j-i). The action of g 
in coordinatized form is then 

(xi, . .. ,x n ) ■ g = (xi ■ 51,... ,x n ■ g n ). 
Killing and Building by Levels 

We call the coordinate tuple the base state, if it consists 
of only the identities (as coset representatives), which 
clearly represents the identity of the original group. 
Given an arbitrary coordinatized state x = (x\, . . . , x n ), 
we call the coordinatewise changes of values from Xi to 
1 (top-down) 'killing by levels'. This is accomplished 
by simply applying the inverse of the coset representa- 
tives, in order. An example is shown in Fig. 01 Con- 
versely, 'building by levels ' is accomplished bottom-up 
by successively applying the coset representatives, i.e. 
elementary generating symmetries whose product is the 
given coset representative Xi for the zth component, in 
the order x n then x n -\ and finally x\ to move from 
the 'solved state' (1, . . . , 1) to create state x bottom-up. 
Moreover, one can compute elements that change only 
a single coordinate to a desired value. 

Global Transformation via Coordinate Values 

Since we work with groups, whenever we make any ac- 
tion, it can be undone by an inverse, thus reversing the 
killing by levels we can go from the base state to any 
other coordinate value combination. Thus going from 
x to y coordinatewise can be achieved simply combin- 
ing the level-killers of x with the level-builders of y. 
More efficient solutions are generally possible (and im- 
plemented), but this provides at least one way to do it 
using the hierarchical coordinate system. 

5 An Application: Permutation Puz- 
zles 

Permutation puzzles are one person games where 
the moves are permutations, elements of a group 
|Joyner, 2002 1. Natural problems for such puzzles are: 



9i+i ■= Xi- gi - [Xi ■ gi 



i+l 



Thus g on x is coordinatized as g(x) = (gi, . . . , g n ). 
Note that generally Xi • gi does not equal to x\ ■ gi, so 



1. How can one go, via elementary legal moves, from 
one configuration x of the puzzle to a standard 
"solved" configuration? 

2. More generally, how can one go from configuration 
x to another arbitrarily selected configuration y? 

The main quest of permutation puzzles is often to find 
a shortest sequence of moves that leads to the solution 



(bounded by the diameter of the Cayley graph of the 
underlying group). We should mention that this is not 
our aim here. We would like to facilitate understanding, 
but usually the quickest solutions are 'dirty tricks' that 
are very difficult to grasp and one has to fall back on 
simple memorizing. However, we suspect that Lagrange 
decompositions might also help in tackling the shortest 
solution problem. Here we focus on 1 . and 2. 

5.1 Coordinatizing Rubik's Cubes 

The 3x3x3 Rubik's Cube is probably the most popular 
permutation puzzle. 

What does it mean 'to know the Rubik's cube'? The 
question usually boils down to the ability to solve the 
cube. By asking a cube-fan he/she would give a few 
tricks, recipes to apply in certain situations. By learning 
these algorithms one can learn to solve the cube, but 
does it imply understanding the cube, i.e. grasping how 
certain sequences of moves work and seeing why they 
work? Not necessarily. 

We claim that understanding comes with imposing a 
coordinate system on the underlying algebraic structure. 
Here we demonstrate this on the symmetry group of the 
Pocket Cube, which is the 2x2x2 version of the Ru- 
bik's Cube. The moves are the 90 degree clockwise ro- 
tations of the 6 sides. By the Frobenius-Lagrange de- 
composition we know that each coordinate system cor- 
responds to a subgroup chain, so devising new algo- 
rithms for solving the cube is equivalent to construct- 
ing subgroup chains. Such a coordinate system encodes 
a 'global viewpoint' in which one solves by successive 
approximation, with manipulations going from coarse 
to fine resolution, and convergencing in terms of mov- 
ing from natural, abstract states to fully specified states. 
The group of the cube acts on the set of configurations 
in such a way that any non-trivial permutation yields a 
different result on the 'solved state', thus the stabilizer 
of this state is (1), so each group element corresponds 
to a unique configuration. Hence by the remark follow- 
ing Theorem [H for any coordinate system arising from 
any subgroup chain down to (1), each configuration of 
the cube has a unique coordinatized form. Examples de- 
rived computationally follow (see also Fig. 01): 

Pocket Cube: Cornerwise Decomposition 

One can solve the cube in a rather long, systematic step- 
by-step fashion: get the position and then the orientation 
of the first corner right, then proceed to the next corner 
until the cube is solved. In the subgroup chain we put 
the stabilizer of the position of a corner, then continue 




Figure 1: A random, 'scrambled' configuration of the 
Pocket Cube, coordinatized by the cornerwise decom- 
position as (8,2,5,3,2,2,5,1,2,3,3,3,1,2). Note: 
coset representatives have been integer-encoded in these 
examples. 
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Figure 2: Pocket Cube configuration after killing the top 
9 levels out of 14. (1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 3, 1, 2) 



with the stabilizer of the orientation of the corner within 
the position stabilizing subgroup. Then we repeat the 
whole process for another corner. The chain yields the 
following coordinatization: 

S s lC 3 lS 7 lC 3 lS 6 lC 3 lS 5 lC 3 lS^C 3 lS 3 lC 3 lC 2 lC 3 

where the top level component S$ acts on 8 
states/coordinate values, representing the 8 possible po- 
sitions of the first stabilized corner. Therefore killing the 
first level will put the corner in the right position. The 
coordinate values on the second level correspond the 3 
possible rotational states of the corner. The 3rd and 4th 
level similarly encode the second corner, and so on. 

Another Model for Understanding: Permute the 
Corners, then 'Beat the Clock' 

Contrasting to the previous, very machine-minded solu- 
tion, here is another one which is short, and reveals the 
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Figure 3: The solved state of the Pocket Cube with co- 
ordinates (1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) 

existence of a different puzzle within the Pocket Cube: 

7 

i=i 

The top level component is the right regular represen- 
tation of the now familiar symmetric group permuting 
the 8 corners. The second level is the direct prod- 
uct of 7 copies of modulo 3 counters (the orientation 
group of corners). It is to be noted that there are not 
8 copies, otherwise every corner could be rotated in- 
dependently from the other corners (and that would 
be rather easy to solve). Actually solving the bot- 
tom level is the same type of problem as the Rubik's 



Clock |Wiggs and Taylor, 1989], which is an array of 
connected modulo 12 counters. As the underlying group 
is commutative, it is easier to solve since the order of 
operations generating this subpuzzle does not matter in 
this lowest level. 

For an example computational session using our de- 
composition package in GAP [GAP, 20061 see Fig- HI 

3x3x3 Rubik's Cube 

Going to the standard Cube we immediately meet some 
difficulty, as its group is not a transitive one. Therefore, 
using (G, G) we can get a decomposition which solves 
the corners as in the Pocket Cube and nearly separately 
and in parallel the remaining non-corner middle faces 
(those not at the corners, not in the middle of a side) on 
which cube group is transitive. Then we can proceed 
by coordinatizing and solving the Pocket Cube and this 
middle cube puzzles independently. 

6 Conclusion and Future Work 



Lagrange coordinatizations of that permutation group, 
as well as to different solving strategies for manipula- 
tion. Coordinatewise manipulation of the permutation 
group via short or minimal length words over group's 
basic generators is an easily achieved next step. For 
exploitation of the idea of hierarchical coordinatiza- 
tion in more general settings, groups can generalized 
to semigroups in order allow the possibility of some 
irreversible manipulations. Computational semigroup 
theory is not yet as advanced as its group theoreti- 
cal counterpart but it is being developed quickly (e.g. 
l Egri-Nagy and Nehaniv, 2008[ ). 
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pocket_cube_F := ( 9 , 10 , 11 , 12 ) ( 4 , 13 , 22 , 7 ) ( 3 , 1 6 , 2 1 , 6 ) ; 
pocket_cube_R := ( 13 , 14 , 15 , 1 6 ) (2 , 2 , 22 , 10 ) ( 3 , 17 , 2 3 , 11 ) ; 
pocket_cube_U := ( 1 , 2 , 3 , 4 ) ( 5 , 17 , 13 , 9 ) ( 6 , 18 , 14 , 10 ) ; 
pocket_cube_L := ( 5 , 6 , 7 , 8 ) ( 1 , 9 , 2 1 , 1 9 ) ( 4 , 12 , 2 4 , 18 ) ; 
pocket_cube_D := (2 1 , 22 , 2 3 , 2 4 ) ( 12 , 1 6 , 2 , 8 ) ( 11 , 15 , 1 9 , 7 ) ; 
pocket„cube_B := ( 17 , 18 , 1 9 , 2 ) ( 1 , 8 , 2 3 , 14 ) (2 , 5 , 2 4 , 15 ) ; 
pocket_cube_gens := [pocket_cube_U, pocket_cube_L, pocket_cube_F , 
pocket_cube_R, pocket_cube_B, pocket_cube_D ] ; 
pocket_cube_gen_names : = [ "U" , "L" , "F " , "R" , "B" , "D " ] ; 
pocket_cube := GroupByGenerators (pocket_cube_gens) ; 

gap> subgroupchain := Chief Series (pocket_cube) ; 

[ <permutation group of size 88179840 with 6 generators>, 

<permutation group of size 44089920 with 10 generators>, 

<permutation group of size 2187 with 7 generators>, Group (() ) ] 
gap> Remove ( subgroupchain, 2 ) ; 

<permutation group of size 44089920 with 10 generators> 

gap> decomposition := LagrangeDecomposition (pocket_cube, subgroupchain); 

Reducing the number of generators 6->2 7->7 DONE 

Lagrange components .. DONE 

Lagrange mappings .. DONE 

Getting component names.. DONE 

Guessing state sets.. DONE 

Generating states. . .DONE 

Coset action mapping DONE 

Lagrange decomposition of : <permutation group of size 88179840 with 6 generators> 

1 40320 <permutation group of size 40320 with 2 generators> 

2 | -2187 <permutation group of size 2187 with 7 generators> 
gap> StructureDescription (decomposition [1] ) ; 

"S8" 

gap> StructureDescription (decomposition [2 ] ) ; 
"C3 x C3 x C3 x C3 x C3 x C3 x C3" 

gap> scrambled := Random (pocket_cube) ; 

(1,21,17,22, 6,3,20,8) (2,16, 9,10,15,19,18,12) (4, 13,23,24,5, 7, 14, 11) 
gap> coordinates := Perm2CascadedState (decomposition, scrambled); 
[ 16074 1438 ] 

gap> levelkillers := LevelKillers (decomposition, coordinates); 
[ (1,24, 15, 3, 4, 11, 2, 12) (5, 19, 20, 13, 9, 16, 17, 21) (6, 22, 14, 7, 18, 8, 23, 10), 
(1, 18, 5) (2, 14, 17) (3, 10, 13) (4, 6, 9) (8, 19, 24) (15, 23, 20) ] 

gap> epi := EpimorphismFromFreeGroup (pocket_cube : names : =pocket_cube_gen_names ) ; 
[ U, L, F, R, B, D ] -> [ (1, 2, 3, 4) (5, 17, 13, 9) (6, 18, 14, 10) , 

(1,9,21,19) (4,12,24,18) (5,6,7,8), (3,16,21,6) (4,13,22,7) (9, 10, 11, 12), 

(2, 20, 22, 10) (3, 17, 23, 11) (13, 14, 15, 16), (1, 8, 23, 14) (2, 5, 24, 15) (17, 18, 19, 20), 

(7, 11, 15, 19) (8, 12, 16, 20) (21, 22, 23, 24) ] 
gap> PrelmagesRepresentative (epi, levelkillers [ 1 ]) ; 
U*B~2*R"-l*F*L*D*L"-l*F"-2*U~2*F*L*F*L"-l*tr-l*F"-l 
*U*F*U"-1*L*F"-1*D*L"-1*F"-1*U*R"-1*U 
gap> PrelmagesRepresentative (epi, levelkillers [2 ]) ; 
F"-1*U*F*U"-1*L*F~-1*D*L~-1*F"-1*U*R~-1*U*F"-1*L~-1*U*L 
*U*F*U"-1*F*L*F*U*L"-1*U"-1*L~-1*U~-1*R~-1*F~-1*L*D*L"-1*B~-1*L~-1 

gap> halfsolved := scrambled * levelkillers [ 1 ] ; 
(1, 5, 18) (2, 17, 14) (3, 13, 10) (4, 9, 6) (8, 24, 19) (15, 20, 23) 

gap> half solvedcoords := Perm2CascadedState (decomposition, halfsolved); 
[ 1 1438 ] 

gap> halfsolved * levelkillers [2 ] = (); 
true 

Figure 4: Deriving a Coordinate System for Pocket Cube and Solving via Killing by Levels. This excerpt from 
a GAP interactive session, using our implemented software package for automatically generating and manipulating 
coordinate systems on permutation groups. The possible moves as generators (basic symmetry operations F, R, U, 
L, D, B rotating a face 90°) for the Pocket Cube are defined together with names. Group elements are represented 
here in cycle notation, with () denoting the identity element (or, according to the right regular representation, the 
solved state). The subgroup chain for the decomposition is the modified chief series of the group. Then the decom- 
position is calculated yielding the two level coordinatization. Then a scrambled (random element of the Pocket Cube 
permutation group) is shown in coordinate format. Finally, the scrambled cube is solved by levels of this hierarchical 
coordinate system, top-down, using level-killers (see text), which are also expressed as an [unoptimized] sequence 
of the original generators. 
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